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ABSTRACT 

Secondary anisotropics of the CMB have the potential to reveal intricate details about the history of 
our universe between the present and recombination epochs. However, because the CMB we observe is 
the projected sum of a multitude of effects, the interpretation of small scale anisotropics by future high 
resolution experiments will be marred by uncertainty and speculation without the handles provided 
by other observations. The recent controversy over the excess small scale anisotropy detected by 
CBI and the BIMA array is a foretaste of potential challenges that will be faced when interpreting 
future experiments. In this paper we show that cross correlating the CMB with an overlapping weak 
lensing survey will isolate the elusive kinetic Sunyaev-Zel'dovich Effect from secondary anisotropics 
generated at higher redshifts. We show that if upcoming high angular resolution CMB experiments, 
like PLANCK/ACT/SPT, cover the same area of sky as current and future weak lensing surveys, 
like CFTHLS/SNAP/LSST, the cross correlation of cosmic shear with the kSZ effect will be detected 
with high signal to noise ratio, increasing the potential science accessible to both sets of surveys. For 
example, if ACT and a CFHTLS like survey were to overlap this cross-correlation would be detected 
with a total signal to noise ratio greater than 220, reaching 1.8 per individual multipole around 
/ 5000. Furthermore, this cross-correlation probes the three point coupling between the underlying 
dark matter and the momentum of the ionized baryons in the densest regions of the universe at 
intermediate redshifts. Similar to the tSZ power spectrum, its strength is extremely sensitive to the 
power spectrum normalization parameter, cts, scaling roughly as (t|. It provides an effective mechanism 
to isolate any component of anisotropy due to patchy reionization and rule out primordial small scale 
anisotropy. 

Subject headings: cosmology: theory ~ cosmology: observation - cosmology: weak lensing - cosmology: 
peculiar velocities - galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe 
(WMAP) 1 has ushered in an era of unprecedented 
accuracy for measurements of the cosmic microwave 
background (CMB) anisotropy. Future high angular 
resolution experiments like PLANCK^, the Atacama 
Cosmology Telescope (ACT)^, and the South Pole Tele- 
scope (SPT)* will measure arc-minute scale (£ > 1000) 
anisotropics at the /iK level. The primary anisotropics 
generated a t the epoch of recom bination and measured 
by WMAP l|Spergel et al. 1 120031 involve calculations in 
linear perturbation t heory an d arc on sound theoretical 
footing (see however iBean et al. I l200ll. and references 
therein for a discussion of alternative models). However, 
secondary anisotropics, caused by highly nonlinear struc- 
tures and involving complicated dissipative baryonic 
physics, are still a subject of theoretical speculation 

In considering secondary anisotropics, a practical 
distinction must be made between anisotropics that 
have a CMB-like thermal spectrum, and the anisot- 
ropics with a non thermal frequency dependence, most 
notably the th ermal Sunyaev-Zel'dovich effect (tSZ) 
ijSunvaev fc Zel'dovich 11980 ). Although challenging and 
imperfect, it should be possible to isolate these ther- 

^ http://map.gsfc.nasa.gov/ 

^ http:/ /astro. estec.esa.nl/SA-general/Projects/Planck 
^ yhttp: / / www.hep.upenn.edu/ ~angelica/act/act.htmli 
* [Tittp: / / ast ro. uchicago.edu/spt| 



mal contributions from a temperature map using the 
specific frequency and spatial dependence of the latte r 
l|Bouchet fc Gisperill99alTegmark fc Efstath iou 11996(1 . 
After these non-thermal components have been removed, 
the dominant sources of secondary anisotropics are ex- 
pected to be — from degr ee to arc-minute scal es — 
the Rees-Sciama effect (R,Sl ()R.ees fc Sciama II1968I) . the 
weak gravitational lensing of the CMB itself, the kinetic 
Sunyaev-Zel'dovich (kSZ) effect and possibly patchy 
reionization (Sunyaev & Zel'dovich 1980 and Santos et 
al. 2003 for recent references). Given a "frequency 
cleaned" temperature map containing the sum of these 
secondary effects, it would be desirable to study each 
component individually in order to evaluate their con- 
sistency with theoretical predictions. If the CMB tem- 
perature were studied in conjunction with another low 
redshift tracer of the cosmological density field, the pres- 
ence or absence of correlations between the two might 
provide a mechanism to isolate the various components. 

The recent detection of excess small scale anisotropy 
is a nice illustration of the preceding discussion and pro- 
vides a foretaste of the potential challenges that will be 
faced when interpreting the secondary anisotropics de- 
tected by future experiments. At angular scales of a sev- 
eral arc-minutes at multipoles beyo nd the dampin g tail 
(^ > 2000), the CBI experiment Mason et al. 1 120031 
and the BIMA array (|Dawson et al. 1 120021) have both 
detected temperature anisotropy at a level of ^ 500/j,K^ 
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with significance of > 3cr. Although a natural interpre- 
tation of this excess power is tSZ, the strength of the 
signal requires that the power spectrum normalization 
parameter take on a value erg « 1 in the upper range 
allowed by current CMB and Large Scale Structure data 
l|Bond eit oi. ]\2W^ iKomatsu fc Seliak 20021. Further- 
more, the required value is likely even higher because 
the measured anisotropy, if thermal SZ, was lik ely di- 
luted by radio point source subtraction iHolder J I 2Q_03 ■ 
Both the CBI and BIMA experiments observe at low fre- 
quencies ^ 30 GHz well into the Rayleigh-Jeans region 
of the thermal SZ frequency spectrum, so at present the 
non-thermal frequency signature of the tSZ effect cannot 
be exploited to determine the nature of this excess. This 
has led to considerable speculation about other sources 
of arc-minute scale anisotropy possibly arising from the 
epoch of reionization or the early universe. While the 
kSZ effect and patchy reionization are not expected to 
produce small scale anisotropics at the level detected, it 
has been suggested that the excess power is due to SZ 
fluctu ations from high-z > 10 star formation (|0h et al. I 
|2003|) . primordial voids create d by bubble nucleation 
durin g the inflationary epoch l)Griffiths. Kunz. fc Silk I 
1200 3') , or broken scale invariance i n the primordial power 
spect rum produced by inflation l)Coorav fc Melchiorri I 
Unn^). Cross correlating the CMB temperature maps 
measured by CBI and BIMA with a low redshift tracer 
of the density field would provide a mechanism to de- 
termine whether the small scale excess was generated in 
the local universe. The detection of a correlation with 
galaxies would favor tSZ or radio source contamination 
rather than the more speculative alternatives. 

Cosmologists will likely face similar uncertainties in the 
future for experiments that will have the frequency cov- 
erage required to produce frequency cleaned temperature 
maps that effectively remove the tSZ. Cross correlation 
with a local density tracer could provide the most ef- 
fective means to determine the source of the anisotropy. 
For example, the patchy reionization signal or some ex- 
otic form of small scale anisotropy generated by infla- 
tion will not correlate with low redshift density tracers; 
whereas, the kinetic SZ effect should correlate strongly, 
which brings us to the subject of this work. Weak grav- 
itational lensing, or the coherent distortion of images of 
faint background galaxies by the foreground matter dis- 
tribution (see Van Waerbcke & Mellier 2003, for a re- 
cent review), is a tracer of the local density field, with the 
added advantage that it probes the dark matter directly, 
foregoing the complications caused by issues of bias in 
galaxy surveys. In this paper we compute the correla- 
tion between weak gravitational lensing and the kinetic 
Sunyaev-Zel'dovich effect and evaluate the prospects of 
future experiments to measure it. 

Although the cross-correlations between secondary an- 
isotropics and weak gravitational lensing of the CMB it- 
self, has been considered previously (Spergel & Gold- 
berg 1999, Goldberg & Spergel 1999, Cooray & Hu 2000, 
Takada & Sugyama 2002, Verde & Spergel 2002), the po- 
tential correlation signals between the CMB and galaxy 
weak lensing, or "cosmic shear" has been somewhat ne- 
glected. In particular, we focus on previously unexplored 
small angular scales beyond the "damping tail," which 
will be probed by the aforementioned future experiments. 

The overlap of angular scales between secondary an- 



isotropics of the CMB and cosmic shear is illustrated in 
figure n Error bars are shown for current and future 
CMB experiments and weak lensing surveys. Theoret- 
ical estimates for primary and secondary anisotropics 
are shown in the left panel, and predictions for the 
power spectrur n of the weak lensing converge nce field is 
shown at right. iSeliak. Burwell. fc PenI l)2001fl calculated 
the correlation between cosmic shear and the thermal 
Sunyaev-Zel'dovich effect (tSZ). Once the tSZ signal is 
removed by frequency cleaning, the Rees-Sciama effect, 
CMB lensing , kin etic SZ, and possibly patchy reioniza- 
tion remain. iHu I |2002) considered the correlation be- 
tween cosmic shear and both the linear ISW effect and 
CMB lensing. As is visible in the figure, the nonlinear 
correction to the ISW effect, known as the Rees-Sciama 
effect is important on small scales £ > 300, and we 
compute the correlation between cosmic sh ear and the 
Rees- Sciama effect in a companion paper l|Dore et al. I 
I2003j) . For £ > 3000, the kinetic SZ spectrum intersects 
the lensed CMB damping tail, and will likely dominate 
the anisotropy spectrum, (we have plotted the maximal 
patchy reionization signal for illustration). 

Below we calculate the fully nonlinear cross correlation 
between cosmic shear and the kinetic SZ effect. This cal- 
culation is complicated by the fact that the kSZ effect 
is proportional to the velocity of the ionized baryons, 
AT/T ~ v; hence, any two point correlation with a den- 
sity tracer will be negligible for an isotropic velocity field. 
This is simply the statement that the over-density field S 
is just as likely to correlate with a cluster moving toward 
us as one moving away from us, and thus the average 
correlation vanishes. In fact, the isotropy of the velocity 
field guarantees that any statistic of odd powers of the 
velocity field is highly suppressed relative to even statis- 
tics (Monin & Yaglom 1971, Scannapieco 2000, Castro 
2003). Thus, we must work with a three point statistic 
- two kSZ points and one weak lensing — and we con- 
sider a "collapsed" configuration. That is we calculate 
the two point function of weak lensing and the filtered 
CMB temperature squared, which condenses three point 
information into an easily measurable power spectrum. 

We model the observables in § |5] and calculate the an- 
gular auto power spectra of weak lensing and kinetic SZ 
which are shown fig. ^ Special attention is paid to the 
approximations and analytical techniques used to cal- 
culate these auto power spectra in the fully nonlinear 
regime. The cross-correlation signals are calculated in 
§ El In § 0] we attempt to qualitatively understand why 
cosmic shear is correlated with the square of kinetic SZ, 
and the dependence of the correlation signal on the power 
spectrum normalization is explored in § El A signal to 
noise ratio analysis is performed in § El to quantify the 
prospect for detecting this correlation for current and 
future experiments. We discuss our results, discuss limi- 
tations, and conclude in § [3 

Throughout this paper wc consider only flat cosmolo- 
gies. As our fiducial mo del we c hoose the best fit WMAP 
(only) model of Sp ergel et al. I Ipl031 with n„i = 0.270, 
r^A = 0.73, h = 0.72, Hs = 0.97, as = 0.84 and t = 0.17. 

2. MODELING THE OBSERVABLES 

We define in this section the needed general equa- 
tions before introducing the effect that will be discussed 
throughout this paper. We calculate their angular auto 
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Fig. 1. — Left: Expected signal and errors for the coming PLANCK and ACT CMB experiments (see JH]for precise assumptions). The 
Bordeaux wine color refers to PLANCK whereas the cream one refers to ACT. The cosmic variance is computed assuming the primordial 
power spectrum only, which is obviously an lower limit. The thick solid line denotes the primordial CMB signal with (up) or without (down) 
the lensing contribution. The thick short-dashed line denotes the ISW effect (low t) and the non-linear RS effect (high i). Note the angular 
transition between both. The thick long-dashed line denotes the expected thermal SZ contribution according to iKomatsu & Seliak 2002). 
The thick double-dot-dashed line corresponding to the expected kSZ, whereas the dot-dashed line corresponds to a possible component 
induced by p atchy-reionisation (note that the overall amplitude is rather model-dependent and that we here assume the strongest signal 
ISantos et al. iji200 3)). The thin solid line denotes the sum of all those contributions. iJig/it'Expected signal and error for the coming 
CFHTLS, SNAP and LSST experiments (see ijS|for precise assumptions). We consider here only shot noise due to intrinsic ellipticities and 
cosmic variance, and assume that both dominates over systematics at any scale. We assume for this illustration only that this 3 surveys 
measure the same convergence power spectrum. We thus neglect the effects of the different sources populations. From darker to lighter, 
the 3 depicted error boxes correspond to CFHTLS, SNAP and LSST. 



power spectra and sketches their cross-correlation. The 
calculations in this section are complicated by the fact 
that the three dimensional fields whose power spectra 
and bispectra we desire are highly nonlinear. Special 
attention is paid to the approximations and analytical 
techniques used to calculate these nonlinear power spec- 
tra. 



2.1. Angular Power Spectra in the Flat Sky 
Approximation 

We will be interested in the angular power spectra of 
fields on the sky X{i\) which are weighted line-of-sight 
projections of three dimensional fields which we denote 
5x, 



X{n) 



where W''^ (rj) is the weight function for the field X and 
f-^ is a real space isotropic filter with which we convolve 
this field. We define the conformal "look back" time as 
r] ^ tq — t, which can also be interpreted as the comoving 
distance. Here r is the conformal time dr = (1/a), tq, 
which can also be interpreted as the comoving distance 
from the observer. Here and throughout we set c = 1 

For small sections of the sky or high multipole mo- 
ments, it is a good approximation to treat the sky as flat. 
In this flat-sky approximation, the Fourier moments of 
this field on the sky are 



X{i) = / d-^n X{n)e 



(2) 



Combining with eq. Q gives 

X{£) = f^{l) J ^W^iv) J ^5(k^=£^,fe)e^'=^^ 

^ (3) 
where f'^{t) is the Fourier transform of f^. 

As in the all sky case, the cross correlation (or auto 
power spectrum) for any two fields, X and Z is defined 

by 

{X*{e)Z{l)) = {2ti)Hd{1 + (4) 

and similarly the power spectrum between two three di- 
mensional fields 5x and 5z is 

(5x(k,77)^z(k,r;)) = (27r)35c(k + k')Pxz(fc,??), (5) 

Where (5_d is the Dirac delta function. 

For the relationship between these flat-sky Fourier 
modes and the spherical harmonic coefficients Xim and 
a proof of the correspondence betwe en the angu lar power 
spectrum Cf ^ and C^ ^ffl, see e.g. IWhite et al I lfl999l) 
or Appendix C ofiH^I l|2000l) . 

Combining eqs. ©, Q, and the correlation be- 
tween X and Z is 



{X*{£)Z{E')) = {2^ff^{l)f{l) / dry 



(rj) 



./ ?7'^ ./ 27r 



Provided that the window functions W{vi) are slowly 
varying, k /W'^ , w e can work in the Limber ap- 

proximation ijKaiser 1 IT991 . allowing us to neglect radial 
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modes relative to perpendicular modes, <^ i/rj, and 
we arrive at 

A|zW = ^/^W/^W I d\nrjTj'w''iv)W^iv) 

xAj,z{k^i/v,v) (7) 
where the we have introduced the dimcnsionless angular 
power spectrum A'j^^{i) — -^Cf^, and the dimension- 
less 3-d power spectrum A'j^^{k,r]) = k^ /27T^Pxz{k,v)- 
In the following section, we list the relevant window 
functions and discuss the phenomenological techniques 
used to calculate the three dimensional power spectra in 
the fully nonlinear regime. Note that until § |S1 we set 

2.2. Weak Gravitational Lensing 

The distortion of background source galaxies by the 
foreground matter distribution is completely described 
by the convergence field ^(n) in the weak lensing ap- 
proximation, which is a weighted projection of the den- 
sity field along the line of sight 



(8) 



K(n) = / d-qW^ir/) S{r]h,'i]), 
Jo 

where the kernel is given by (for a flat cosmology) 



2a 



dr]' 



,v_n 



Sill') 



(9) 



and S{if) describes the normalized radial distribution of 
sources 

Sir,)=pAz)^, (10) 
az 

with Pz{z) being the normalized source redshift distribu- 
tion J dzpziz) = 1, which in principle can be measured 
from the weak lensing survey. For definiteness, we use 



Pz (z) 



1 

24 



(11) 



This redshift distribution peaks at 2zo, has mean red- 
shift (z) = 3zq, and has been used in previous studies 
of cosmic shear (Wittman et al. 2000). The values of 
zo that we use to represent current and future surveys, 
along with the other specifications are listed in table 2. 

In order to evaluate the weak lensing power spec- 
trum we need an expression for the fully nonlinear 
power spectrum of the over-density field Pssik). Fit- 
ting formula for the n onlinear power spectruni have been 
studied extensively llHamilton et aT] I199H iMa 1 119981: 
IPeacock fc Dod ds "199( iHSmith e t al. ^20m. We use the 
scaling formula of Smit h et al. I (2003) for the nonlinear 
power spectrum, which is accurate at better than the 
10% level. 

The angular power spectrum Cg of the cosmic shear 
is shown in fig. H Error bars are shown for current and 
future surveys which will have different source redshift 
distributions, although for illustration, we show only the 
case Zo = 0.5. 

2.3. The kinetic SZ effect 

CMB photons diffusing to us from the surface of 
last scattering are up-scattered by the electron plasma 



bound in the gravitational potentials of dark matter ha- 
los. The motion and density variations of those scat- 
terers will imprint a specific thermal temperature fluc- 
tuation pattern on the CMB, O""^^ = ATusz/Tq, both 
because of the Dopp lcr effect and of variations in opacity 
^Sunvaev fc Zel'dovich 1980). 



e'^S^(n)= / dr? g (77) n • p(fi77, 77) 

Jo 



(12) 



where we have introduced the momentum density 

p(firy,7y) = (1 + (5(n77,?7))v(n77,77) (13) 

and the visibility function, g. The latter can be expressed 
in terms of the Thomson optical depth, t[vi) — T{rf)drj, 
where f{rf) = a(ri)neXeO'T with ctt the Thomson scat- 
tering cross section, Ue the electron density and the 
ionization fraction. Thus 
dr 

g{v) = -r^-^ = ^eTH{l + zfe-\ (14) 



where 



drj 



TH = 0.0691 (1 - Yp) Qbh, 



(15) 

is the Thomson scattering optical depth to the Hubble 
distance today, assuming full hydrogen ionization, and 
Yp is the primordial helium fraction. We assume further 
that reionization occurs instantaneously, i.e. = 1 at 
all times after the reionization epoc h, Zri. Then we have 
iCxriffiths. Barbosa. fc Liddig[T99l 



r[Zri) — —— — 
3 i^jji 



(16) 



for a flat universe. 

As was first pointed out bv iKaiserl l|1984^ . cancella- 
tion of successive peaks and troughs for small wavelength 
modes prevents a significant kinetic SZ effect from con- 
tributing to the CMB anisotropy at linear orde r. This 
can b e understood in the context of the 'Lim ber' ijKaiser I 
IT99I or 'weak coupling' iHu & Whitel[T99l approxima- 
tions. They dictate that only modes with wave-vectors 
perpendicular to the line of sight contribute to the pro- 
jection of a given random field, because uncorrelated ra- 
dial waves cancel in the line of sight projection. Gravity 
generates potential velocity flows, so that radial veloc- 
ities have no wave-vector component perpendicular, in 
Fourier space, to the line of sight; consequently, the lin- 
e ar Doppler effect p ^ v van ishes. 

lOstriker fc Vishniad l)1986D showed that the same will 
be not be true for density modulated velocity flows, 
p ^v, which can have a non- vanishing curl. Thus 
wave-vectors perpendicular to the line of sight gener- 
ate fluctuations parallel to the line of sight. At sec- 
ond order in perturbation t heory, this is known as 
the Ostriker-Vishniac effect llPodelson fc JubasI Il995t 
iHu fc White! 119961 l()st'riker fc Vishniad 119861: iVishniac I 
1987); whereas, the full nonlinear signal is referred to 
as the kinetic-SZ (kSZ) effect. This full non-linear 
contri bution has been stud ied in d etail bot h analyti- 
cally llHull200nHMafc Frvl2001: Zha.ng ^^^,112003) and 
numerically fda Sil va et 0/71120011 iGnedin fc Jaffai200ll 
[Hpringcl et al. 20011: iZhansf et al 1120031) " 

Although the amplitude of the expected signal is un- 
certain, the maximal signal can be calculated if we as- 
sume that the electrons trace the dark matter exactly. 
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Then we can write for the 
momentum component pa 



power spe ctrum of the radial 
Psl (17) 



p±p 



(fc) « -I'l-nis 



where ^'pj^p^ is the power spectrum of the vortical com- 
ponent of the momentum field p^. This assumes the 
dominant contribution to the power spectrum of the ra- 
dial momentum comes from nonlinear densities coupling 
to linear bulk velocity flows (Hu 2000; Ma & Fry 2001; 
though see also Zhang et al. 2003). The right panel 
of fig. 01 shows the 3-d auto power spectrum Pp-p- at 
z — 0.5. The projection of this power spectrum with 
eq. gives the angular kSZ power spectrum, which is 
shown in fig. H 

3. THE KINETIC SZ WEAK LENSING CORRELATION 

In this section we compute the correlation between the 
kinetic SZ effect and weak gravitational lensing. This cal- 
culation is complicated by the fact that the simple two 
point correlation will vanish because Q^^"^ ^ v, 

and the velocity field is isotropic. This is just the state- 
ment that the over-density field 5 is just as likely to cor- 
relate with a cluster moving toward us as one moving 
away from us, and the average correlation vanishes. The 
isotropy of the velocity field guarantees that statistics 
of odd powers of the velocity field are highly suppressed 
relative to even statistics (Monin & Yaglom 1971, Scan- 
napieco 2000, Castro 2003). 

Thus a non- vanishing correlation requires that we work 
with a three point statistic: two kSZ points, so that our 
statistic is even in the velocity, and one point weak lens- 
ing. We consider the simplest case of a "collapsed" three 
point function — the cross correlation between the kappa 
field and the square of the kSZ, — which condenses 

three point information into an easily measurable angu- 
lar power spectrum. 

Squaring the temperature field couples multipoles at 
all scales, and in particular, power from multipoles out- 
side the range where the kSZ is dominant, will swamp 
the signal. To prevent our quadratic temperature statis- 
tic from being polluted by the beam and primary anisot- 
ropics, we must filter the temperature before squaring 



(18) 



This simple convolution before squaring is a special case 
of a more general class o f filt e rs for quadratic tempera- 
ture statistics studied bv lHu 1 1)2002() . 

3.1. Limber approximation for the bispectrum 

In the flat sky approximation, the square of the filtered 
temperature fluctuation can be written as a convolution 
in (2D) Fourier space. 



eK^') 



{24 



(19) 



and the weak lensing kinetic SZ squared correlation is 



^ I ±l^{ki£)Qfinefi£'-t)) (20) 



with 



rj^ 



rj'- 



(21) 

27r 7 27r ^ ' 

X Bsp,p,{£/v + k„r/r/' + k^', {£' - £")/»/ - (k, + k^')), 
where we have introduced the hybrid bispectrum 

CS (ki)pfl (k2)pfl (k3)} = 

(27r)3 Bsp.p, (ki, k2, kg) Sd (ki + k2 + ka) . (22) 

The expression in eq. H22I) can be simpl ified in the Lim- 
ber a pproximation for the bispectrum (^Buchalter et aLl 
|2000|) . valid for the small angles considered here. We can 
again ignore radial modes kz ^ (/rj, giving 

( mOfinQfii' - i")) « i27r)^Sn{£ + ^')/(^")/(l^' - ^"D 

>^ / [5 Bsp.p, {Ih. tM, {£' £")/?;) (23) 

Plugging this into eq. H20|l gives 

{K{i)e}{i')) = {2TTfSD {i + e') c^Hi) 



(24) 



with 



drj 



C ' W = j ('?) [9 {V)] nk = i/v, (25) 

and where we have defined 

'^ik,ri)= I — —2 /((777)/(|k+q|77)B5p.p. (k,q, -k- q) , 
J (27r) 

(26) 

This power spectrum in eq. H25() can be thought of 
as a Limber projection of the "power spectrum" T(fc,7y) 
which is an integral over all planar triangle conflgura- 
tions of the bispectrum (k, q, — k — q), with one of the 
sides of length k. It depends only on the norm of k by sta- 
tistical isotropy and effectively collapses the three point 
information contained in S^p.p. (|Coorav Il2001bj) . We 
henceforth refer to eq. (|26|l as the triangle power spec- 
trum. 

3.2. The hybrid bispectrum Bsp-pf^ 

In order to calculate the weak lensing kinetic SZ 
squared correlation, we must evaluate the hybrid bispec- 
trum i?5pfipfl of the density and radial momentum com- 
ponent. Here we proceed by analogy with previous an- 
alytical studies of the kSZ power spectrum. The power 
spectrum of the radial momentum component is 

(Pfl (k) pn (k')) = {27rf 5d (k + k') Pp^p, (fc) . (27) 

iHu I ll2000[l and then the joint numerical and analytical 
work of iMa fc Frvl l)2000D suggested that the dominant 
contribution comes from non-linear densities coupling to 
linear bulk velocity flows. Furthermore, because of the 
cancellations that occur in the projection (Kaiser 1984), 
only the curl component of the projected momentum p 
- i.e. the component perpendicular to n in the Fourier 
domain that we thus note p_L- will contribute to the kSZ 
temperature fluctuation, so that 



{mQfinQfi^'-n) = {2nrmf{\£' -£"\) 



rPr. 



(fc) 



7,2 pnl 



(28) 
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where v'^^^ is the volume averaged velocity dispersion 



(29) 



We follow this analytical line of reasoning and make 
the following ansatz 

1 o „„, ^3Q^ 



B, 



Spapa ~ q ^nns^SSS i 



which presumes that the dominant contribution to the 
bispectrum in eq. (|22|l comes from large scale bulk 
velocity flows coupling to the three point function of 
the nonlinear density field. Note that at the 2 points 
level, analogous ap proximations can be moti vated in a 
halo model context ijCoorav fc Sheth 2002; S heth et al. I 
1200,'^ . Note however that it was also claimed, even if 
not clearly illustrated, that the curl component of the 
momentum induced by the vorticity in the non-linear 
velocity field itself might not be completely negligible 
ijZhang et al. \ t2003). If this were of some importance, 
then our evaluation would underestimate the signal of 
interest. 

We are interested in the high-A;, non-linear behav- 
ior of T{k), thus we need to evaluate the integral in 
eq. (|26|l with the fully nonlinear bispectrum of the den- 
si ty field. We use for this purpose the fitting function 
ofi coccimarro &: Couchmani l|200H) for the bispectrum, 
i.e. 

Bf,, = 2Ff (ki, k2) (fci) ih) + cyclic (31) 
where 

(ki,k2) = ^a{ncs,ki)a{ncff,k2) 



1 ki • ka fk 



2 kik2 ki 



b{ncs,ki)b{ncs,k2) 



2 /ki -kj 



c(ricff,fci)c(neff,A;2) • (32) 



7 \ fcifc2 

The fitting functions are given by 

1 + cTs-°-2 (z) x/0.7 Q3 M (g/4)"-+^-^ 



a (neff, k) = 
6(noff,/c) = 

c(noff,/c): 



1 + ((7/4)""«+3-' 
1 + 0.4 (noff -I- 3) g"<=«+3 

4.5 \ /'o„^"cff+3 



(33) 



1 + (2g) 



^^cff +3.5 



where rics is the effective spectral index of the power 
spectrum, 

,,, dlnP 
^cff {k) = . (34) 

The quantities q, Q, and as (z) are defined by 

with |ipiin(fc„,,z) = l, (35) 



4 _ 2"'=tt 
1 + 2""ff+-'- ' 



(36) 



and 



<78 (z) 



(2^ 



P]^"(fc,z)W^(fci?) (37) 



where W{kR) is the usual Fourier transform of a top-hat 
of radius R = 8 h^^ Mpc. 

The functions a(ncff,fc), b{nes,k), and c{ncs,k) in 
eq. (|34ll interpolate between the one loop perturba- 
tive expansion and highly non-linear regimes for general 
CDM cosmological models. It can be seen that for large 
scales, i.e. k <^ k^i, a = b = c — 1 and the tree level 
perturbation theory expression is recovered, whereas on 
small scales, i.e. k ^ fcni, a = cr^"'^ (z) Q3 (ricft) 

and 6 = c = 0, so that the bispectrum becomes indepen- 
dent of triangle configuration. 

3.3. The Triangle Power Spectrum 

In this section we evaluate the triangle power spec- 
trum for the simplest case where the filter functions are 
set to unity, f{£) = 1. In this case, a simple approxima- 
tion exists in the high k limit that dramatically simplifies 
the computation and provides and intuitive understand- 
ing of the source and strength of the correlation. The 
full computation will be performed in for the purpose 
of computing the signal to noise ration in SJS] one the 
appropriate filter is introduced. 

The integral in eq. H26() will have three terms corre- 
sponding to the three permutations of the wave vectors 
in eq. 

r (k) = [Ti (fc) + T2 (fe) + % (fc)] (38) 



where 



Ti (fc) = 2P (fc) I Ff (k, q) (q) (39) 

T2(fc) = 2P(fc) / -f^Ffik,-k^<DP^U\^ + <l\) 
J (27r) 

%ik)^2 [ ^F2°«(q,-k-q)P,t(g)P,t(|k + q|) . 

We first remark that T2 — Ti after a simple transla- 
tion of the integration variable. Then, 7i can be eval- 
uated exactly in polar coordinates because the angular 
integration factorizes out of the radial integration 



77r 



5a{k) I dqqa (q) Pj'i (q) 



c{k) / dqqc{q)Pf^{q) 



(40) 



where we abbreviate a {n^s, k) as a (fc), and likewise for 
b and c. The remaining term 7^, can be written 



%{k) 



d\ 
(2^ 



a{q)ai\k + q\) 



fe(q)fe(|k + q|) 
c(q)c(|k + q|) 



q| q 

2 



^)(^") 



p,"i(g)p,"i(|k- 



-q|) 

(41) 



with /i = k ■ q. 

As mentioned previously, we arc interested in the high- 
k, non-linear behavior of T (fc) as the weak lensing kinetic 
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SZ squared correlation will only be significant on small 
angular scales. The dominant contributions to the inte- 
grals in eq. (|4U|) will be for q near the peak of P^l {q) 
where a[q) ^ b (q) ^ c{q) ^ I . Thus in the high-/c 
limit the integrand is significant only for q k, and 
we can drop terms of order O (q/k). A similar approx- 
imation has been used previously to calculate the mode 
coupling integrals for the kinetic SZ power spectrum (Hu 
2000a; Ma & Fry 2001; Cooray 2001). In this approxima- 
tion 73 = 7i, and we get, introducing the dimensionless 
power spectra, and Aj-^j 

k^ 

A^- = — r and Al 
the following simple result 



^nl 



h.3 

^ pnl 



Al (fc, z) = -vl^^ (z) A^i (fc, z) E3 (fc, z) 



(42) 



(43) 



where we defined the three point enhancement of the 
power spectrum 



E-- 



5a{k,z) / d\nq a{q,z) 



A^i {q,z) 



+ c{k,z) / d\nqc{q,z) 



(44) 



Here the time dependence has been explicitly included 
in the functions above as a reminder. With these defini- 
tions, we obtain the final expression for the weak-lensing 
kSZ squared power spectrum 

^^.e^ ^^Jdr,r^W^ (,7) [g {rj)f A^- (fc = i/r,, r,) . 

(45) 

This intuitive result states that is given by 

a 'Limber' projection of the triangle power spectrum, 
A|-(A;,z), which condenses information from all planar 
triangle configurations of the hybrid density-momentum 
bispectrum with side length k. The form of the triangle 
power spectrum in eq. (|43|l has as simple interpretation. 
It is similar to the kSZ power spectrum in eq. H28I) . in 
that it arises from density modulations of a large scale 
bulk flow. However, here we are dealing with three point 
modulations, which are enhanced by the non-linear cou- 
pling of small scale modes to large scale power. This 
mode coupling is encapsulated in the integrals over the 
power spectrum in E-^ {k, z). The dimensionless triangle 
power spectrum eq. (|43ll is shown in the right panel of 
fig. along with the kSZ auto power spectrum eq. H28|l . 
Note that they have similar orders of magnitude. 

4. UNDERSTANDING THE CORRELATION 

Before calculating the signal to noise for the aforemen- 
tioned cross correlations, we attempt to understand the 
degree to which the secondary anisotropies are correlated 
with weak lensing. 

The strength of the angular correlation between any 
two random fields X and Y is quantified with the cross- 
correlation coefficient 



Corr-^'*'(^) 



XY 



YY 



(46) 



As we are dealing with angular correlations of pro- 
jected three dimensional fields (see eq. JT))), this cross- 
correlation coefficient will depend both on the extent to 



which the window functions VF^(ry) and {i]) overlap, 
and the strength of of the cross power spectrum Pxy 
relative to the auto power spectra Pxx and Pyy- 

If we change the integration variables from rj to k — 
£/rj, we can rewrite eq. Q for the angular power spec- 
trum, as an integral over wavenumber 



A' 



where 



XY{i)= dink Hi{k)Aj,z{k,V^i/k) 



Hnk)^^i^^j W''{£/k)W''{£/k) (48) 

With the above two equations, we can think of any cross 
correlation as a weighted integral of the cross power spec- 
trum, A^j-y, with the projection 'kernel' Hf^{k). The 
kernels Hf^{k) and Hj^ (k) indicate the scales probed 
by X and Y respectively, while the cross kernel Hf^{k) 
indicates their degree of overlap. 

The cross correlation coefficient for the kSZ^-weak 
lensing correlation is plotted in fig. 13 Computing this 
required the auto power spectrum of the kSZ tempera- 
ture fluctuation squared C^^^" (see §5 eq.( [5^ 1. The 
level of correlation is significant, approaching ^ 0.8 at 
the arc-min angular scales £ > 3000 in the damping tail, 
where the primary anisotropies are heavily attenuated. 

The physical scales probed by these power spectra is 
illustrated in flg. 01 where we plot the projection kernels 

He{k), for C;«, Cf^'^^^, and for £ = 1000 and 

10000. The 3-d power spectrum of the radial momentum 
Ap.p. and the triangle power spectrum A^ are plotted in 
the right panel at 2 = 0.5. Note that the triangle power 
spectrum is of the same order of magnitude as the power 
spectrum responsible for the kSZ effect, A^ ~ ^pnPa- 

The angular power spectrum is the projection of the 
power spectra on the right weighted by the kernels on 
the left (see eq.(EIl), in this flgure. The area under the 
Hi{k) are normalized to unity. The sharp cutoff in the 
kernel for Cf''^-^^'^ , occurs because we have assumed 
the universe is reionized abruptly at rj{z-cci — 20) so that 
W^^'^{£/k) drops to zero at small k. 

We show the differential redshift contribution to each 
angular power spectrum as a function of redshift for £ — 



(47) 



dz 



for the auto 



5000 in flg.0 where we plot 

and cross spectra. By deflnition, the area under each 
of these curves integrates to unity. While the kSZ auto 
power spectrum (dashed) has contributions from a high 
redshift tail, the bulk of the signal is skewed toward low z 
~ 0.5, because of the growth of nonlinear structure. The 
overlap with the weak lensing (dotted) is thus apparent, 
explaining the large correlation coefficient found in fig. El 
(left). 

5. DEPENDENCE ON POWER SPECTRUM NORMALIZATION 

It is well known that the thermal SZ power spectrum 
exhibits a strong dependenc e on the power spectrum 
normalization parameter cts ("Komatsu & Seliak ' l2002t 
IHpliak. BurwcU. & Pen 2001; Zhang & Pen 2001), ap- 
proximately scaling as the seventh power Ci oi a\. This 
strong dependence suggests that the thermal SZ power 
spectrum might be an effective probe of cxg , because even 
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Fig. 2. — Le/J.Angular power spectra for kinetic SZ and kSZ^-weak lensing cross correlation. Only dimensionlcss units are used. The 
blue solid line corresponds to the kSZ^ auto power spectrum, while the red solid line corresponds to the WL auto power spectrum scaled 
down by a factor 10^ for the sake of illustration. The kSZ^-weak lensing angular cross power spectrum is shown by the solid (black) line. 
The dotted black line shows the result of using the approximation in eqn. I4iji for the triangle power spectrum, which is good to ~ 25% 
for i > 3000. The dashed (blue) line is the kinetic SZ auto power spectrum Right: Cross correlation coefficient for the kSZ^-weak lensing 
correlation. 




Fig. 3. — Projection Kernels for the weak lensing kinetic-SZ correlation: The left panel shows the projection kernels Hi (fc) for the the 
weak lensing auto correlation (dotted) C^", the kinetic-SZ auto correlation (dashed) C^^^"**^^, and the weak lensing-kinetic SZ squared 

cross-correlation (solid) C^^"^ All kernels have been arbitrarily normalized to unit area. The leftmost set of curves indicate the range 
of wave-numbers contributing to the angular power spectrum at £ = 1000; the rightmos t se t is for i = 10000. The right panel shows 
dimensionlcss 3-d power spectra of the radial momentum component (dashed) . (eq. I28i and the triangular power spectrum (solid) 

(eq. 1431 . The triangle power spectrum is of the same order of magnitude as the power spectrum responsible for the kSZ effect. 



if theoretical predictions are off by a factor of 2, this 
translates into less than 10% systematic error in erg. 

We can determine the rough erg dependence of the 
kSZ^-weak lensing correlation from simple power count- 
ing if we recall that C'^^z -k^ ^ v'^^^Bsss- The matter 
power spectrum scales as Pgg oc tTg"^ — two powers of 
(jg in the linear regime and three in the strongly nonlin- 
ear stable clustering regime. The large scale modes that 
give rise to the bulk flows in v'^^^ are linear, so we ex- 
pect w^ms PsT '^8- Finally, the bispectrum scales as 
Bggg cx Pgg , so putting everything together we expect the 

scahng C^^^ oc Cg"®. Figure [S] illustrates the depen- 
dence of the cross power spectrum on erg. We find that 



^ksz -K ^ ^7 pi-ovides a good description of this scal- 
ing. This is similar to the strong seventh power scaling of 
the thermal SZ and suggests that the kSZ^-weak lensing 
correlation might also serve as an effective probe of the 
power spectrum normalization and its redshift evolution. 

6. PROSPECT FOR MEASUREMENT: SIGNAL TO NOISE 
ANALYSIS 

In order to assess the detectability of the kSZ^-weak 
lensing correlation we consider several current and future 
CMB experiment configurations and WL lensing survey 
specifications 

On the CMB side, we consider the future PLANCK 
mission and the soon to be build ACT ijKosowskv Il2003(l 
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Fig. 4. — Differential redshift contributions to angular power 
spectra: Differential redshift contributions to the weak lensing 

(dotted), the kinetic-SZ auto correlation 

(dashed), and the weak lensing-kinetic SZ squared cross- 

" at £ = 5000. All curves have been divided by 
the total power so that they integrate to unity. 



auto correlation ^C: 
correlation t^C', 

Z7r ' 




100 



1000 

I (multipole) 



10000 



Fig. 5. — Dependence of the kSZ^ — k cross correlation on erg. 
From top to bottom, the lines indicate ag=1.2, 1.1, 1.05, 1.0, 0.95, 
0.9, and 0.8 as labeled in the figure. 



telescope. Whereas the first one should offer a full sky 
survey (/sky — 0.85 if one takes into account galactic con- 
tamination) the second one should conduct a smaller (100 
square degree (/sky — 0.0025)) but deeper survey. Both 
precise characteristics are summarized in table 1, where 
^fwhm denotes the beam full width at half-maximum, api^ 
denotes the instrumental noise per 0'^ area pixels and 
"Area" denotes the observed sky area. 

The subsequent expected errors are illustrated on fig.H 
Both have a frequency coverage appropriate for a proper 
tSZ "separation" so that we neglect this signal from now 
on (but the effects of residuals in § . On the lensing 
side, we consider 3 various surveys: the on- going Canada 
France Hawaii Telescope Legacy Survey (CFHTLS) ^, 

^ |http://www.cfht.hawaii.edu/Science/CFHTLS| 



the future SNAP satellite ^ and the future ground based 
Large Synoptic Survey Telescope (LSST) . Their key 
characteristics are summarized in table 2 and the ex- 
pected uncertainties are illustrated in fig. ^ : zq de- 
notes the redshift parameter of the source distribution 
of eq. rigai the mean density of galaxies, the sur- 

veyed area and the single component rms shear due to 
intrinsic ellipticity is a-y. 

We define the optimal signal to noise per individual 
mode £ as 



Cov[(C^ 



XY\ 



(49) 



To derive this formula, the assumption was made that 
all the relevant fields were Gaussian so that we end up 
naturally with the S/ {S + N) Wiener optimal weighting. 
Doing so, we follow e.g. Zaldarriaga ( 2000). We then 
approximate the covariance matrix by its diagonal, i.e. 
neglecting the correlation induced by the cut sky, 



(/sky(2^+l)) 'C0V[(C/ 

{cf+Nn{cr+Nr) + {cry 



(50) 



where Nf^ denotes the evaluated instrumental noise 
contribution to the measured angular power spectrum 
of £. For the weak lensing, it is simply the shot noise due 
to the intrinsic ellipticitie of the sources, i.e. 



(51) 



whereas for the CMB, we include in the noise definition 
all the signal that can not be separated out using multi- 
frequency information and the conjunction o f the noise 
and beam smearing give rise to l|Knox Ill995() 



I instr 



fwhm 



(52) 



In the following subsections, we will consider two 
types of overlapping surveys. Either a small survey, like 
CFHTLS cross-correlated with either ACT or SPT, in 
which case the overlap area is limited by the sky cover- 
age of say ACT, which is 100 square degrees, or a very 
large survey which uses PLANCK in conjunction with 
LSST or PANSTAR. In the latter case the total over- 
lapping sky is limited by LSST, which is around 30,000 
square degrees. 

6.1. Signal to noise for the kS^-WL correlation 



Table 1 

Main characteristics of the CMB experiments 

DISCUSSED. 





^fwhm 




Area 






[mK] 


fdeg.^l 


ACT/SPT'' 


1.7 


2. 


100. 


PLANCK 


5. 


2.2 


34000. 



"SPT has in fact a wider but shallower survey but the correlation 
would thus be limited by lensing survey. 



http://sn ap.lbl.gov| 
nttpT/^™™ • l^^t o . org! 
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1 10 100 1000 10000 1 10 100 1000 10000 



Fig. 6. — Signal to noise estimation for logarithmic bins in I: The left panel assumes an ideal full sky measurement free of any instrumental 
noise. The middle panel illustrates the expectation for a 100 square degree of a CFHTLS like survey and ACT. The right panel assumes a 
LSST like 30 000 degree square degree and PLANCK. For each of those plot, the dashed line corresponds to a worst case scenario where 
the patchy reionisation contribution is 10 times higher (in temperature square units) than the kSZ one. 



For this correlation, the dominant "noise" contribution 
that need to be considered at the relevant angular scales 
(~ 2', i.e. I ~ 5000, see fig.^l are the lensed primordial 
CMB, the potential patchy reionisation contribution and 
the instrumental noise. We neglect any residual tSZ sig- 
nal and the subsequent correlation it would induced, we 
also neglect the "spurious" correlation due to lensing of 
the primordial CMB, although we discuss both of them 
in the next section. As is evident from fig. ^ the am- 
plitude of spurious correlation with lensing of the CMB 
should be be more than 2 orders of magnitude weaker. 

To evaluate both C\^'^ "^^^ and involves one 

more approximation. We indeed have to evaluate the 
power spectrum of quantities squared in real space. To 
do so exactly would involve the computation of even 
higher order statistics, e.g. it would be a 8 point func- 
tion for C^^^ . To avoid such complications, we 
can reasonably consider those fields as being weakly non- 
gaussian, so that we can neglect this way their connected 
part, in particular the trispectrum of the squared field 
|Zaldarriaga 2000). With this hypothesis, the evolution 
of the power spectrum of those quantities translates in 
a convolution in harmonic space that can be easily per- 
formed numerically : 

C^\C)^2 J ^^C'ii') C^{t~t') (53) 



Table 2 

Main characteristics of the weak-lensing surveys 
discussed. 





zO 


"gal 

[gal.arcmin^^] 


Area 
[deg.2] 




CFHTLS 


0.5 


20. 


170. 


0.31 


SNAP 


0.9 


100. 


300. 


0.23 


LSST 


0.75 


75. 


30000. 


0.16 



^ j C^(|.g-£'|) 

We evaluate the convolution above numerically using 
the full power spectrum of primary plus secondary an- 
isotropics. With this, we can compute all the required 
quantities in eq. (|49|l . 

Before computing the theoretical prospects for the ex- 
perimental configurations we consider, it is crucial to re- 
alize that before squaring the Q field, an appropriate fil- 
tering is necessary. Indeed, as visible e.g. in eq. (|53|l . 
this squaring will introduce violent mode-coupling, so 
that any signal that has a strong power outside the range 
where the kSZ signal dominates will pollute the quadratic 
field. It would be possible using our theoretical knowl- 
edge of the shape of the cross-correlation power spec- 
trum C'^^^^^ , although a bit cum bersome, to compute 
an optimal filter for our statistics l)Hu 1 120021. Instead 
we take a reasonable short cut. Given the fact that both 
CI and C^^''^^ are nearly flat in the range where the 
kSZ dominates, we will consider indeed as the filter for 
our cross-correlation signal the one that maximizes the 
kSZ signal alone. This last one can be easily derived if 
we assume that the projected kSZ signal is weakly non- 
Gaussian. It is indeed given by the usual Wiener filter 
l|Bouchet fc GisDerdlTooa: iTesfmark fc Efstathiou 1119961) 

r'ksz 

^ CNoisc (7kSZ ' 

where C^°"''^ includes the instrumental noise but also all 
the other astrophysical component different from the kSZ 
itself. Note that this filter is similar to the one derived 
in a more general manner by Hu (2002) and employed 
bv lCoorav l l )2001a|) . 

The results are illustrated in fig. El where we plot the 
contribution to the square of signal to noise per log £. We 
consider two situations: i/ the dashed line corresponds 
to a situation where the patchy reionisation is present 
an d strong, i. e. we coii sider the higher amplitude model 
of ijSantos et al. Il2b03|) . where the patchy reionisation 
power spectrum amplitude is about 10 times higher than 
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with weak lensing at the two point level. This would 
effectively "lensing clean" the temperature map so that 
the residuals would be dominated by the kSZ. Realistic 
simulations are certainly required to precisely evaluate 
the potential of this technique, but this is outside the 
scope of this paper. 

7. CONCLUSION 

We have evaluated the prospects for upcoming exper- 
iments to measure the correlation between arc-minute 
scale secondary anisotropics of the CMB and cosmic 
shear, after tSZ has been removed by frequency clean- 
ing. In particular, we evaluated the signal to noise ratio 
for the kSZ^-weak lensing correlation. 

The two point correlation of kSZ with and density 
tracer is negligible because of the isotropy of the velocity 
field; however, we found that a strong 3-point correlation 
exists. A collapsed three point statistic was introduced to 
measure the three point signal, which is the angular cross 
power spectrum between weak lensing and the square of 
the filtered temperature. While residual thermal SZ will 
limit the detectability of this correlation for PLANCK- 
LSST, it should be easily detectable by cross correlation 
ACT/SPT with the CFHTLS/SNAP at ^ - 5000 where 
the signal to noise ratio per ^ reaches 1.8. In principle, a 
full analysis of the cross bispectrum B^^^^^ could yield 
a higher signal to noise — a calculation we defer to a 
future paper. 

The kSZ^-weak lensing correlation probes three point 
correlations between the underlying dark matter and the 
the momentum of the ionized baryons in dense regions, 
providing the only known probe into this physics. We 
presumed that these higher order correlations, or more 
precisely the hybrid bispectrum Bsp^^p^^ , arose from three 
point density modulations of large scale coherent mo- 
tions. The situation is surely more complicated and this 
approximation must be compared to hydrodynamic nu- 
merical simulations. Nevertheless, measurement of kSZ^- 
weak lensing correlation will provide valuable insights 
into the physics of energy injection, the bias and ion- 
ization fraction of baryons in the densest environments, 
and information on the nonlinear mode coupling between 
dark matter and baryons. We also considered the de- 
pendence of correlation signal on the power spectrum 
normalization parameter as and determined the approx- 
imate scaling cx cr| similar to the scaling of the tSZ 
power spectrum. 

Despite the interest in detection of the kSZ^-weak lens- 
ing correlation, a null detection could prove equally in- 
teresting. The contribution of patchy reionization to arc- 
minute scale anisotropy in the CMB is a contentious 
issue, with different reionization models and histories 
predicting sig nal amplitudes differing by several orders 
of magnitude JAehanim. Desert. PuEet. & GisDertlllQQd 



the kSZ one in the same angular range (see fig. ii/ 
the solid line corresponds to a scenario where the patchy 
reionisation component is completely negligible. 

Obviously, this signal should be easily detectable us- 
ing a CFHTLS-ACT configuration. The overall signal to 
noise ratio reaches 225 if the patchy reionization is neg- 
ligible, 47 otherwise. The signal to noise per individual £ 
reaches 1.8. For the PLANCK-LSST configuration, the 
signal reaches 978 without no patchy reionization and 
473 otherwise. However, as we discuss in the next sec- 
tion, this kSZ signal will probably be heavily contami- 
nated by thermal SZ residuals. 

6.2. Spurious Correlations 

The primary limitation to detecting the kSZ^-weak 
lensing correlation will certainly come from spurious 
sources of correlation. Possible sources of anisotropy 
other than kSZ that the weak lensing may correlate 
with are tSZ residuals from imperfect frequency cleaning, 
weak lensing of the primordial CMB, and point sources. 
We discuss each in turn. 

The dominant spurious correlation will probably be 
that due to t h e resi dual tSZ after frequency cleaning. 
ICoorav et aTl ll200(]() estimated that for PLANCK, the 
power spectrum of the residual for 600 < £ < 2500 should 
be around 0.14 times the initial tSZ power spectrum, 
i.e. around 5 times the kSZ for our model. Since the 
correlation coefficient around WL and tSZ is exp ected 
to be around 0.6 l|Seliak. Burwell. fc Penll2001|) . the 
correlated component between tSZ residuals and WL 
will probably be more significant than the correlation 
with kSZ. For ACT/SPT telescope, the frequency cov- 
erage will allow for frequency cleaning as good as if not 
better than PLANCK, while the higher resolution pro- 
vides improved sensitivity at the scales where the tSZ 
dominates, £ > 2000 (see fig. QJ, so frequency cleaning 
should be drastically improved at those scales. Further 
consideration of this issue is required, but we defer this 
to a future work. 

Another source of spurious correlation should be the 
correlation induced between the lensed primordial CMB 
and the WL. Here again, this effect might be important 
for PLANCK as visible in fig. ^ but should be negligible 
for ACT/SPT since in the range of interest, the lensed 
contribution of the CMB is around 2 orders of magni- 
tude weaker than the kSZ one. In addition, any att empt 
to cr oss correlate CMB lensing with cosmic shear ()Hu I 
120021) will need to consider the kSZ^-weak lensing signal 
as a source of spurious correlation, since this signal will 
begin to dominate for £ 3000 near the damping tail. 

Finally, there is the possibility of spurious correlation 
with point sources. At the frequencies of interest, this 
is primarily dust emission from sub-mm galaxies. This 
population of objects will probably correlate with the 
weak lensing signal but the strength of that correlation 
is rather uncertain. 

Finally, we note that the kSZ^-weak lensing correlation 
has a specific signature that may alleviate some of the 
problems with spurious correlation and could provide a 
mechanism to isolate the kSZ signal. Specifically, the 
kSZ effect does not correlate with weak lensing at the 
2 point level, whereas tSZ residuals or dusty sub-mm 
galaxies will. It may be possible to exploit this fact, 
and remove the part of the anisotropy that correlates 



'Gruzinov & Hu ' '1998'; 'Knox. Scoccimarro. fc Do delsoiJ 

1998; Santos et al. 2003; Valagcas et al. 2001). The 
lack of a correlation in the presence of significant arc- 
minute scale anisotropy would indicate that this anisot- 
ropy was generated at high redshift, rather than in the 
low redshift universe probed by weak gravitational lens- 
ing. This last point is especially important when one 
considers the similarity in shape between the power spec- 
trum of patchy reionization and that of kinetic SZ (see 
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In addition, we note that many of the calculations in 
this paper can be directly applied to correlations with 
galaxies rather than weak lensing as the tracer of den- 
sity field, given a suitable model for the bias of galaxies 
in the nonlinear regi me, as provided b y, e.g. t he recently 
popu lar halo model ijCoorav fc Sheth 2002; S heth et al. I 
12003(1 . The weak lensing window function W'^{vi), would 
simply be replaced with a suitable window function for 
the galaxies. As the noise in the correlation measure- 
ment is dominated by the CMB, the signal to noise 
per square degree should not change significantly. How- 
ever, because large area galaxy surveys already exist, 
the correlations we discussed might be measured with 
gala xies first, as was the case recently for the ISW ef- 
fect ( Afshordi et al. "2003; 'Boughn & Crittenden "200j 
iFosalba fc E G atzanaga 2 003; Fosalba et a l. 20dl 
INolta et al. l T2003: Scranto n"^oi. 1 120031) . The kSZ^- 
galaxy correlation would probe three point correlations 
between the density of galaxies and the momentum of 
ionized baryons in clusters, providing valuable informa- 
tion about galaxy formation. 

With such a high signal to noise ratio, it is conceivable 
that radial information might also be extracted from the 
cross-correlation. If photometric redshifts of background 
source galaxies were available, one could attempt to de- 
project the cross-correlation tomographically. Similarly, 
if galaxies were correlated with kSZ^ rather than weak 
lensing, galaxy photo-z's could be use to extract this ra- 



dial information. This has been d iscussed in the context 
of galaxy-tSZ cross corr elation bv lZhang fc Pen 1 12001). 
Studying the cross-correlation as a function of redshift 
would yield information on the redshift evolution of the 
dynamical state of baryons in dense environments. 

In conclusion, our study predicts that coming arc- 
minute scale measurements of the secondary anisot- 
ropics of the CMB should will correlate strongly with 
local tracers of the density field, like weak gravitational 
lensing or galaxies from large scale structure surveys. 
The total signal to noise ratio for the kSZ^-weak lens- 
ing correlation is greater than 220 for ACT/SPT cor- 
related with CFHTLS/SNAP, which will be easily de- 
tectable. We believe that this signal along with the 
other cross-correlation s with cos mic sh ear studied by 
ISebak. Burwell. fc PenI ll20(Tl and IhTTI {2002), provide 
a significant incentive for upcoming CMB experiments 
and lensing surveys to image the same regions of sky. 
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